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o 

The linear dispersion relation of a finite amplitude, parallel, circularly polarized Alfven wave in a 
relativistic electron-positron plasma is derived. In the nonrelativistic regime, the dispersion relation 
has two branches, one electromagnetic wave, with a low frequency cutoff at yl-f- 2lj| / 0,% (where 
Lu p = (4irne 2 /m) 1 ^ 2 is the electron/positron plasma frequency), and an Alfven wave, with high 
frequency cutoff at the positron gyrofrequency £l p . There is only one forward propagating mode 
for a given frequency. However, due to relativistic effects, there is no low frequency cutoff for the 
electromagnetic branch, and there appears a critical wave number above which the Alfven wave 
ceases to exist. This critical wave number is given by ck c /Q, p = a/77, where a = wH/f2p and 77 is the 
ratio between the Alfven wave magnetic field amplitude and the background magnetic field. In this 
' case, for each frequency in the Alfven branch, two additional forward propagating modes exist with 

Oh, equal frequency. 

A simple numerical example is studied: by numerically solving the coupled system of fluid and 
Maxwell equations, normal incidence of a finite amplitude Alfven wave on an interface between two 
electron-positron plasmas of different densities is considered. 

! I. INTRODUCTION 

o 
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Electron-positron plasmas are different from electron-ion plasmas, because in the absence of a mass difference, 
there are no high or low natural frequency scalesjll Such plasmas are found in pulsar magnetospheres,0 models of 
primitive Universe, active galactic nuclei jets, [J and laboratory and tokamak plasmas. 0,0] Relativistic effects 
are expected to play an important role in several of these systems. Understanding interactions between waves and 
relativistic electron-positron plasmas is relevant to proposed pulsar emission mechanisms 0, , and may give insight 
into structure formation in the early Universe 0. 

Therefore, wave propagation in relativistic electron-positron plasmas has been the subject of many studies, either in 
, the fluid or the kinetic treatments: linear waves [ 1 Oi . 1 1 1| . nonlinear waves HIES, and nonlinear decays El 13 E3 • 
' In this article we deal with an Alfven wave propagating along a constant magnetic field in a pair plasma. When 
fully relativistic effects are considered in the particle motion, the dispersion relation exhibits unique features which, to 
t-H ' our knowledge, have not been discussed before. We then outline the numerical strategies we are currently considering 
\ to examine the consequences of such features. 

O i II. DISPERSION RELATION 

• 1— I 

>> 



- 1—1 

X 



We assume that the electron-positron plasma is described by the following set of equations: 



5h " V dt J I x J J ' m j \ c 



-^ = -V-(n i iT i ) , (1) 
3 + Vj ■ (7^) = — (e + -v, x B) , (2) 



V • E = 4vrp , (3) 
Vx^=-if, (4) 
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Vx B -^/ + If, (5) 
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J = Yl q 3 n ^3 > ( 6 ) 

P = Yl q i n i ' ( 7 ) 
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where nj is the density of each fluid, Vj is the bulk velocity of each fluid, E and B are the electric and magnetic fields, 
respectively, J is the total current, m is the particle mass, and c is the speed of light, j = p for positrons, and j = e 
for electrons. 

We assume that a circularly polarized Alfven wave propagates along the z-axis, as well as the existence of a constant 
magnetic field in the same direction, Bqz. The wave electric and magnetic fields are given by: 

B = B[x cos(fcz — Lot) + y sin(fcz — Lot)] , (9) 



E = —B[x siu(kz — cot) — y cosffcz — Lot)] . (10) 
ck 

Introducing JHJ) and (|10|l in the fluid equations, the transverse velocity for each species is obtained as [lij 

Vn LO Q.A , 

where flj = qjBQ/rrijC is the cyclotron frequency, and rj = B/Bq. The dispersion relation for the Alfven wave is 

r 2 h 2 LO 2 

^T = l-E { 3 (12) 

LO 2 ^> Uip/jU - ilj) 

where loj — (ATrnojqj I ' m j) 1 ^ 2 is the plasma frequency and ncy is the rest density of species j. 

We want to numerically solve the dispersion relation, which is equivalent to simultaneously solving the set of 
equations (JHJ , 1111(1 and H12[l . For an electron positron plasma, lo p =u> e . We define 

lo ck tol v 

x= oT' y= W a= o2' u = ( 3) 

O bp J bp J L- 

Normalized equations (JHJ) , l|ll|l and (|12[1 , are 

7i = (l-^)- 1/a , (14) 
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2 2 a2 ' aa; 

2/ = a: — r , (16) 

X7 p - 1 X7 e + 1 

where the upper (lower) sign in (|15|l is for positrons (electrons). 
Eliminating Uj, the following equation is obtained for 7: 

7? + ^(^) + 7?(-l + i-^) + 7 J ( ± |)-i=0. (17, 

For a given y, from (|17fl 7j is calculated as a function of x, and then the right hand size of Eq. l|16f) can be plotted. 
This is shown in Fig. Since y is given, we plot the left side of Eq. (|16l) on the same graph, and the intersection 
points of the curves are the roots of the dispersion relation. The curve on the right corresponds to the light wave 
branch of the dispersion relation, and the curve on the left corresponds to the Alfven branch. It is interesting to note 
that when y > y c — a/rj the Alfven wave ceases to exist. 

Thus dispersion relation (|16|l can be solved. The result is shown in Fig. El F° r V = [Fig. Ufa)], the nonrelativistic 
result is recovered. There are two branches, one corresponds to the light wave, with a cutoff at lo/Q p = \/X + 2a. The 
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FIG. 1: Graphical representation of 1161 . for a = 1, y = ck/Q, p — 1, rj = 0.2. Full circles: right hand size of 116H : dashed line: 
left hand size of 11611 



other one is the Alfven branch, which has a resonance at the positron gyrofrequency. For rj = 0.1 [Fig.[3Jb)], however, 
the light wave branch has no low frequency cutoff, and the Alfven branch stops at y c = a/rj. For even higher values 
of the wave amplitude rj [see Fig.^c) for r) = 1], the Alfven branch exists for a very short wavenumber range. The 
Alfven branch is also constrained to a shorter frequency range. This can also be seen in Fig. From Eqs. (|14|> — (|lfc>|> 
it can be shown that the Alfven branch in Fig. ^ has a frequency cutoff at the critical frequency 
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FIG. 2: Dispersion relation 11611 for (a) a = 1, rj = 0; (b) a = 1, rj = 0.1; (c) a = 1, rj = 1. Since the relevant parameter is 
y c = a/rj, Figs, (b) and (c) are reobtained if rj — 0.1 is constant, whereas a is decreased from a = 1 to a — 0.1. 



III. OUTLINE OF THE NUMERICAL STUDY 

Fig-IHshows that, for a given frequency w, the existence of Alfven waves and the number of wave modes depend on 
the physical parameters of the plasma. This poses a number of new questions. For instance, we could consider the 
problem of normal incidence of a plane wave on a density interface. In the usual problem of a light wave in a dielectric 
medium, only one mode exists for a given frequency. Given only the boundary conditions that the fields must satisfy, 
it is possible to completely solve the problem and obtain explicit solutions for the field amplitude at both sides of the 
interface. However, now several modes exist for a given frequency. Moreover, if an Alfven wave exists with a given 
frequency in one side of the interface, then if the density is small enough on the other side it could not exist. Does it 
become an evanescent wave? Is it converted to other modes? 

In order to investigate this we are studying the system by means of a fluid simulation, and by direct time integration 
of the evolution equations. 
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For the fluid simulation, we notice that the system equations lead to the following wave equations for the electro- 
magnetic field: 



The wave equations for the fields l|19|) and (|2U|) , and the momentum equation J2Jl , are discretized using time-centered 
and space-centered finite differences. 

The second approach we are developing is time integration of the evolution equations CEJ-© by means of the 
rationalized Runge-Kutta method. 

Both procedures are under development/testing. A particularly simple problem is that of a pure Alfven mode 
propagating in a relativistic plasma. Since in a pair plasma there is no harmonic generation, all oscillatory fields and 
velocities are purely transverse. We are currently working to obtain satisfactory results for this problem, in order to 
later consider incidence on a normal interface, where longitudinal oscillations may not be neglected. 



Propagation of a finite amplitude Alfven wave in an electron-positron plasma has been studied. Full relativistic 
effects on the particle velocities in the wave field have been considered. The dispersion relation for propagation 
along a constant magnetic field has been obtained and numerically solved. Several features are different to the usual 
nonrelativistic result [Fig.^a)]. In the nonrelativistic case, there are two branches, an Alfven wave and a light wave. 
The light wave has a low frequency cutoff at the frequency u>/fl p — \/l + 2a, and the Alfven wave has a high frequency 
cutoff at the positron gyrofrequency, u> — fl p . However, in the relativistic case the light wave branch does not have a 
low frequency cutoff. As to the other branch, there are two Alfven wave modes for any given frequency, and there is 
a high wavenumber cutoff at ck c /£l p = a/rj. Besides, the high frequency cutoff for the Alfven wave depends on the 
plasma parameters. For large amplitude wave or small plasma density, the Alfven branch is constrained to a small 
region of frequencies and wavenumbers [see Fig. EJc)] . 

We intend to investigate the behavior of the system when a relativistic Alfven wave of frequency uj propagates 
through a pair plasma, and then finds a region where plasma parameters are such that such wave cannot propagate. 
The existence of more normal modes than in the usual nonrelativistic problem make this a nontrivial matter. As 
a first example, we are considering the simple problem of normal incidence of an Alfven wave on a planar density 
discontinuity. Two approaches are being developed, namely a fluid simulation, and time integration of the evolution 
equations. Work on this is in progress, and we expect to present results elsewhere. 

It is worth noting that for waves in the high-wavenumber end of the Alfven branch in Fig. particles have large 
Lorentz factors, and therefore synchrotron radiation emission will be important. Thus, the results presented here may 
be further modified when additional effects such like synchrotron loss are taken into account. 
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